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FOREWORD 

The initial effort under this study contract entitled "Lubrication 
Effect on Bearing Thermal Performances" was performed by Dr. Howard 
Yen. Without his assistance, and that of Mr. Jim Chiou, this docu- 
ment could not have been produced. 

Although the program described herein was in its current stage of 
development in June 1977, funds were not available to complete the 
documentation and completely check out the program. Subsequently, 
a superior program was released by SKF Industries , Inc . , for ana- 
lyzing bearing systems. Therefore, this contractual study was 
concluded under continuation funding by using the SKF program en- 
titled "SHABERTH" to determine its applicability to cryogenically 
lubricated bearing systems. 
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INTRODUCTION 


The "Ball Bearing Heat Analysis Program" (BABHAP) is an attempt 
to assemble a series of equations, some of which are non-linear 
algebraic systems, in a logical order, which when solved, provide 
a complex analysis of load distribution among the balls, ball 
velocities, heat generation resulting from friction, applied load, 
and ball spinning, minimum lubricant film thickness, and many 
additional characteristics of ball bearing systems. These equa- 
tions come primarily from Reference 1, but the film thickness 
originated from References 2 and 3. Although high speed digital 
computers have been utilized to simulate certain phenomena in 
high speed ball bearings, it is deemed that BABHAP is the first 
attempt to solve for all critical performance parameters for high 
speed ball bearings within a single computer program. Therefore, 
on verification by test data, BABHAP will provide a unique tool 
for high speed ball bearing design applications. 

Although initial design requirements for BABHAP were dictated 
by the core limitations of the PDP 11/45 computer (approximately 
8K of real words with limited number of instructions, the program 
dimensions can easily be expanded for large core computers such as 
the UNIVAC 1108. The PDP version of BABHAP is also operational 
on the UNIVAC system with the exception that the PDP uses 029 
punch and the UNIVAC uses 026. A conversion program was written 
to allow transfer between machines. 

A complete listing of BABHAP is contained in the Appendix. 
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TECHNICAL DISCUSSION 

BABHAP is designed to simulate the performance of a series of 
single row, angular contact ball bearings mounted on a shaft, 
with either the inner or outer raceway rotating. The analyst 
using BABHAP should be familiar with the properties of bearings 
and lubricants, and have some knowledge of solution techniques 
for non-linear algebraic equations. A copy of Reference 1 would 
also enhance the user's understanding of this document and of 
BABHAP. A limited introduction to the solution of non-linear 
algebraic equations is contained in the PROGRAM BRIEF section 
of this document, describing the function of the NLINEQ routine. 

The analyst must supply several input parameters which describe 
the physical properties of the bearing system and the lubricant. 
Additional inputs are the loads which the bearing system carries. 
The total list of inputs is described in the INPUT PARAMETERS 
section . 

Figure 1 contains the major sections of BABHAP. As with any set 
of non-linear algebraic equations, an initial guess is necessary 
for their solution. Equations (1) , (2) , and (3) provide the 
initial guess for" the relative axial displacement, the relative 
radial displacement, and the misalignment angle respectively. 
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Equations (4), (5), (6), and (7) are intermediate calculations, 
f. = ^i f do 


(5) 
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Figure 2 contains a typical ball bearing configuration describing 
the value of ball position angle iff, which is determined within the 
program based on the number of balls in the bearing. Figure 3 
describes ball center positions which are related to equations 
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With the currently determined information, equations (10) through 
C17) , which form a non-linear algebraic system, are now solved by 
the Newton-Raphson or Method of Steepest Descent technique, for 
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Where X j is 1 for outer race control and 0 for inner race control, 
0A is determined from equation (18) or (19), (Wp/^) . from equation 
(20) and (Wmy^) . is determined from Table 1. ^ 
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Where ^T' is defined by equation (21) . 

tT' = D/a„ (21) 
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Table 1. (Wm/^^) . versus Raceway Control and Rotation 
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Empirical data 1^(6 shown that the balls will roll on one raceway 


and roll and spin on the other. The raceway on which pure rolling 


only occurs is called the controlled raceway. Equations (22) and 


(23) supply the load deflection factors. 
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where equations (24) , (25), (26), and (27) are intermediate calcu- 
lations . 
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Equations (28) and (29) must be solved to determine the auxiliary 
pareuneters a*, b*, and S * from Table 2. 
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0 
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14.839 
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0.9989156 

17.974 
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0.2106 

0.9994785 

23.55 

0.16442 

0.17167 

0.9998527 

1 

37.38 

0.13050 

0 

0.11995 

0 


TABLE 2 

a*, b* and S* as a Function of F(f°) 


A Lagrangian interpolation is performed due to the non-linearity 
of the table. 


It can be seen that although there are but 8 non-linear equations 
to be solved simultaneously, equations (10) through (29) must be 
solved simultaneously. This creates perturbations in the solution 
technique and does not always provide a perfect solution for each 
ball. However, the results are generally deemed adequate. 
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The preceding calculations have been necessary to establish data 
for calculation of the relative axial and radial displacements 
and the misalignment angle for each bearing. The force and moment 
equilibrium is also necessary for this solution. Equations (30) 
through (34) must be solved simultaneously to provide this informa- 
tion. 
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where, again, Xj is 1 for outer raceway control and Xj is 0 for 
inner raceway control. And [ or | | means the values in 

[ ] or I I are evaluated at bearing number h. is the total 

number of balls in bearing h, and equation (35) defines Xj^j . 

Xij = (fi - 0.5)D + S ij (35) 

For a shaft/bearing system of NTS supports shown in Figure 4 , 
equations (30) through (34) represent a system of 5 x NTS equa- 
tions with c^ah' “^rh' ^h» ^rh» ^h h = 1, 2, ... NTS, 

unlcnowns which are solved simultaneously. 


If the solutions to equations (30) through (34) for Sc( r S and 
6 are not sufficiently close to the original guess, equations 
(8) through (35) are solved again until the old and new values 
agree . 

The normal ball load is calculated by equations (36) and (37). 
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The ball contact areas are related by equations (38) through (41). 
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A check on the validity of the assumed race control is performed by 

evaluating equations (42) through (44) . 

E = a*^ (42) 

2S^ 

s = a/b (43) 

Qm Am Em cos (^^i - c<q) > Qn E^ (44) 

for outer raceway control , m = 0 , n = i 
for inner raceway control , m = i, n = o 

The a, b, and Q are for the most heavily loaded ball in the bearing. 
If the inequality is not true, then the assumed race control is 
switched and the solution begins again with equation (1) . If the 
test is not met for either race control, the problem is ill condi- 
tioned due to the input data and the program terminates. If equation 
(44) is valid, then the cage speed, ball spin to roll.cnatioeaiidhh^ll 
rotational speed are calculated by equations (45) through (47) , 

respectively. 
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The ball normal stress is calculated from equations (48) and (49) . 


r: 

U ij TTa.. 


aijbij 


■ ■ ifijf 


1/2 


( 48 ) 



= 

^aojboj 



(49) 


where x and y are zero at the geometrical center of the ball. The 
shaft deflection calculation from equation (50) is for only one ' 
load, one torque and one moment arm. The value for x is given by 
equation (51) . 
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The bearing stiffness is given by equation (52) . 
= 4.77 X 106 Zh D^/2 cos^/^^ ^rh^^^ 
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The minimum lubricant film thickness is related by equation (53) . 
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where c and >7 are given in Table 3, and equations (54), (55), and 
(56) are intermediate calculations. 
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0.065 

0 .548 
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0.088 

0.620 

0.5 
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TABLE 3 

Lubricant Coefficients 


The film thickness is affected by heating and is related by equa- 
tions (57), (58) and (59). 
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The internal heat generation within the bearing can now be solved. 

The viscous friction torque is calculated from equation (60) or . (61). 
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Single Row 


Mist. Lub. 
1 


Grease Lub. 
2 


Jet Lub, 
4 


TABLE 4 

Lubricant Torque Factor 

The torque due to the applied load is determined from equations 
(62) and (63). 
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For angular contact ball bearings, 

p s: 0.00115 
q « 0 . 33 

and F 0 is the larger of the following two equations. 

F yj = 0 . 9 F^j^ cot o( ~ 0 . 1 Fj-j,j 
= Frh 

The torque resulting from ball spinning is related by equations 
(68) through (73). 
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The resulting heat generation from friction, applied load and ball 
spinning is given by equations (74) through (76), respectively, with 
the total heat given by equation (77) . 
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Although the purpose of BABHAP was to determine internal heat genera- 
tion, a virtual dictionary of ball bearing behavior is described by 
the program. 
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NOMENCLATURE 


A concerted effort was made to use unique mathematical symbols and 
abbreviations. However, the complexity and magnitude of the scope 
of these equations led to some duplication. Usually the duplication 
will be for an insignificant intermediate calculation. 
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a 




di 

do 

dm 

D 



E 

fi 

fo 

Fah 

Frh 

fl 

Fs 


Semi -major axis of the most heavily loaded 
ball of a bearing 

Semi -major axis of j^ ball inner raceway contact 

Semi-major axis of j^ ball outer raceway contact 

Table 2 
figure 4 

Semi -minor axis of the most heavily loaded ball 
of a bearing 

Semi -minor axis of j^h ball inner raceway contact 

Semi -minor axis of j^h ball outer raceway contact 


Table 2 

Basic static capacity 

Viscous heating correction factor, equation (58) 
Inner raceway diameter of a bearing 
Outer raceway diameter of a bearing 
Bearing pitch diameter, equation (6) 

Ball diameter 

Modulus of elasticity of ball 
Modulus of elasticity of bearing 
Modulus of elasticity of shaft h 

fl - 1 -^2^ I -1 

I + IFE i r J ' modulus 

Equation (4) 

Equation (4) 

Axial reaction force of bearing h 
Eadial reaction force of bearing h 
Equation (63) 

Equation (66) or (67) 

Static equivalent load 


Units 

in 

in 

in 

in 

in 

in 

in 

lb 

in 

in 

in 

in 

psi 

psi 

psi 

psi 

lb 

lb 

lb 

lb 
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Units 


FCj 

BD 

F(/0)ij 


F(i°) 

hr. 


03 


Wh 


Hyo 

«sj 

^total 

J 

k 

Koj 


Kij 

Jh 

L 

in 

Mv, 


M 




M, 


9D 


M 


y 

Msj 

n 

^roll 


Centrifugal force of j^ball 

Distance between raceway groove curvature radii 
Inner Raceway curvature difference of ball 
Outer raceway curvature difference of ;fh ball 
Table 4 . 

Minimum film thickness without the consideration 
of viscous heating 

Minimum film thickness 

Heat generation due to applied load 

Heat generation due to friction 

Heat generation due to ball spinning 

Total heat generation 

Moment of inertia of shaft h 

Mass moment of inertia of ball 

Thermal conductivity of lubricant 

Load-def lection factor, equation (22) 

Load-deflection factor, equation (23) 

Span of shaft No. h. Figure 4 

Load factor, equation (59) 

Mass of ball 

Reaction moment of bearing h 
Torque due to applied load 
Gyroscopic moment of ball 
Viscous friction torque 
Torque due to ball spinning 
rpm of rotating raceway 
rpm of the cage 
Ball rolling speed 
Equation (64) 


lb 

in 


in 

in 

BTU/hr 

BTU/hr 

BTU/hr 

BTU/hr 

in4 

in-lb-sec^ 

BTU/hr/ft2*> 

Lb-in^*5 

Lb-in^'^ 

in 

Ibm 

in-lb 

in-lb 

in-lb 

in-lb 

in-lb 

rpm 

rpm 

rpm 
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.1.. 

Units 

Ph 

Concentrate load to the shaft. Figure 4 

lb. 

q 

Equation (65) 

— 

Q 

Normal load of the most heavily loaded ball in 
a bearing 

lb 

Qij 

Normal load between inner raceway and ball 

lb 

Qoj 

Normal load between outer raceway and ball 

in 

Ri 

Radius of the locus of inner raceway groove 
curvature center 

in 


/ equation (54) 

in 

s 

a/b 

— 

Sij 

aij/bij 

— 

Soj 

aoj/^oj 

— 

®xj 

Equation (8) , Figure 3 

in 

j 

Equation (9) , Figure 3 

in 


Applied torque. Figure 4 

in - lb 

u 

1/2 (Uball ^raceway^ 

in/sec 

Uball 

Ball linear velocity 

in/sec 

^raceway 

Raceway linear velocity 

in/sec 

Vj 

Equation (10) ~ equation (17), Figure 3 

in 

Wj 

Equation (10) ~ equation (17) , Figure 3 

in 


(fi - 0.5)D + 5'ij 

in 

Zh 

Number of balls in bearing h 

— 


Initial contact angle prior to loading. Figure 3 

rad 

d^o 

Free contact angle 

rad 


Temperature coefficient of viscosity 

1/°R 


Pressure viscosity component 

rad 

<^ij 

• th 

Contact angle between xnner raceway and j ball , 

Figure 3 

rad 

oCoj 

Contact angle between outer raceway and j ball. 

Figure 3 

rad 
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Units 


ball attitude angle 

rad 

y. 

D/dm 

— 

Yij 

Equation (26) 

— 

iToj 

Equation (27) 

— 

Lh 

Bearing stiffness, Equation (52) 

lb/ ft 


Relative axial displacement 

in 


Relative radial displacement 

in 


Dimensionless contact deformation. Table 2 

^ — 


Contact deformation between ball and inner 

raceway 

in 

Soj 

Contact deformation between ball and outer 

raceway 

in 


Dimensionless contact deformation at inner raceway 


Dimensionless contact deformation at outer raceway 

^ah 

Relative axial displacement at bearing h 

in 

^rh 

Relative radial displacement at bearing h 

in 

u 

Coefficient of friction (sliding) 

— 

f^o 

Lubricant viscosity at ambient temperature 

Ibm/f t • sec 

yo 

Kinematic viscosity of lubricant at ambient 
temperature 

ft^/sec 

h 

Poisson's ratio of ball 

— ■ 

'^2 

Poisson's ratio of raceway 


e 

Misalignment angle 

rad 


Misalignment angle of bearing h 

rad 

w 

Angular speed of rotating raceway 

rad/ sec 

Wm 

Angular speed of cage or orbital speed of balls 

rad/ sec 

Wr 

Ball rotational speed 

rad/ sec 

Wroll 

Angular speed of ball rolling 

rad/sec 

Ws 

Angular speed of spinning motion 

rad/sec 

{T^)ij 

+- v» 

Inner raceway curvature sum of j ball 

18 
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xn 



Units 


^ij 

Eq j 


Outer raceway curvature sum of ball 
Position angle of ball 

Elliptic integral of 2nd kind, equation (70) 
Elliptic integral of 2nd kind, equation (71) 

Parameter 


in 


-1 


rad/sec 
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Figure 1 . 


Major Sections of the Ball Bearing Heat Analysis Program 
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Figure 2. Typical Ball Bearing Position Definition 
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curvalure center Cxcd 


Figure 3. Positions of Ball Center and Raceway Groove Curvature Centers 
at Angular Position ui, with and without Applied Load 



Figure 4. Multiple Bearing-Shaft System 
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PROGRAM BRIEFS 


This section consists of short concise descriptions of the function 
of each routine in BABHAP. 



ROUTINE MAIN 


This routine (overlay) is the director for BABHAP, and calls 
all sub-programs (overlays on the PDF 11/45) . Some equation 
evaluations are made in the routine to reduce the number of 
overlays necessary to fit on the PDP machine. All input and 
calculated data are transferred from MAIN to each overlay through 
labeled COMMON. The following diagram describes the solution 
sequence for BABHAP. 



SOLUTION SEQUENCE FOR BABHAP 













SUBROUTINE INPUT 


This routine (overlay) reads all input data and converts it to 
the proper units. A page header is printed, followed by the in- 
put parameters. The routine is structured to reduce repetitive 
input for additional cases and to allow for chosen parameters 
to be set in the BLOCK DATA routine if desired. 

SUBROUTINE INIT 

Initial values are established for the non-linear equations in 
routine (overlay) INIT. The logic for specifying which raceway 
is controlled is also a part of INIT. 

SUBROUTINE BLK4 

The BLK4 routine (overlay) solves for the position of the inner 
raceway groove curvature center after the effects of the applied 
load, relative to the fixed outer raceway groove curvature center. 


SUBROUTINE BLK5 

The supportive equations necessary to find the relative axial 
displacement, relative radial displacement, and the misalignment 
angle are solved in the BLK5 routine (overlay) . 

SUBROUTINE BLK7 

The calculations for normal ball load and the contact area for 
each ball comprise the BLK 7 routine (overlay) . 
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SUBROUTINE BLK9 


The cage and ball rotational speeds are calculated in routine 
(overlay) BLK9. The ball normal stress due to the inner and 
outer raceway point of contact are also calculated. iBearing 
stiffness is also determined in BLK9. 

SUBROUTINE BLKIO 

The primary function of routine (overlay) BLKIO is to calculate 
the lubricant film thickness for the most heavily loaded ball. 

The film thickness is also evaluated for the inclusion of viscous 
heating on the lubricant. 

SUBROUTINE BL1112 

The BL1112 routine (overlay) calculates the viscous friction 
torque, torque due to the applied load, and torque as a result 
of ball spinning. The torque terms are then utilized to calcu- 
late the heat generation within each ball bearing, as a result 
of friction, load and ball spinning. 

BLOCK DATA 

The BLOCK DATA routine (overlay) establishes table values for 

l^bel COMMON variables. The routine can be used to initialize 

y 

or provide default input values for the more common materials, 
reducing the user input list. The structure of the input routine 
will allow overriding default data without program modifications. 


28 


i 



SUBROUTINE NLINEQ 


This subroutine (overlay) solves non-linear algebraic equations 
of the form 

f(x) =0 (1) 


by combining the techniques of the Newton Raphson method with the 
method of steepest descents The derivatives are approximated by 
differencing and the elements of the Jacobian are determined at . 
each iteration 


^k<X) (2) 

A correction vector is determined by solving the linear system 
of equations 

n 


£ 

j=l 


*^kj §j 


-fj,(X) 


k=l, 2 


.n 


( 3 ) 


Then X is replaced with (X +X?) where X initially is a guess, 
and X is determined by a search process to make the estimate 
(X+ X better than the estimate X. 

The criterion for success of the search is 

I(X +X S) < I(X) 

If the program predicts a singularity in the path of the solution, 
the method of steepest descent is used. This method allows the 
convergence to a stationary point which may not be a solution? if 
this occurs, another initial guess should be tried. 

As with all Newton or steepest descent algorithms, the step size 
and initial guess are critical in finding solutions efficiently. 

•The step length, -A, is revised internally to be large or small to 
obtain a solution. However, the step length is not allowed to get 
smaller than the specified data input DSTEP. The following diagram 
describes the basic program flow. 
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Logic Diagram for NLINEQ 
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PARAMETER LIST 


N 

X 

DSTEP 

DMAX 

ACC 

MAXFUN 

IPRINT 

F 

AJINV 

W 

JC 


the number of simultaneous equations to be solved 

One dimensional array of the unknowns 

step length minimum bound 

maximum step length bound 

specified accuracy for residuals 

limit for number of iterations 

"0" for answer only, “I" for each iteration 

one dimensional array of function values 

two dimensional array for inverse Jacobian 

. 2 

one dimensional working space 2N + 5N 
one dimensional working space 
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SUBROUTINE CALFUN 


This subroutine defines the non-linear set of equations. 

F{X) = 0 

to be solved by the NLINEQ routine. CALFUN resides in the BABHAP 
library on the PDP 11/45. The parameters N, X, F are the same as 
in subroutine NLINEQ. 

SUBROUTINE GJR 

This is a supportive routine for NLINEQ, and resides in the BABHAP 
library on the PDP 11/45. GJR solves simultaneous linear equations 
by using Gauss- Jordan elimination with column pivoting. The param- 
eters are supplied by NLINEQ and therefore deleted from this docu- 
mentation. 


SUBROUTINE LAGRAN 

This subroutine is supportive and resides in the BABHAP library on 
the PDP 11/45. LAGRAN uses the lagrangian interpolation formula 
to find the value of 

Y = F(X) 

The technique uses a skip search to accelerate finding the desired 
data, and provides accurate results for non-linear tables. 

INPUT PARAMETERS 

The input parameters for BABHAP have primarily been limited to the 
physical properties of the system - shafts balls, raceway, and 
lubricant. The remainder of the parameters are related to the 
applied loads. Although there are 26 possible input cards, BABHAP 
has been coded to allow the analyst to select as few as one input 
card to run repetitive cases, or to specify any or all inputs in 
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the BLOCK DATA routine without modifying any other segments of 
the program. Figure 1 is a sample list of the input data. The 
required formats and variable names are given in Figure 2. 





Figure 1. Sample Run Stream 
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Figure 1 . Sample Run 


Stream and Data Input for BABHAP (Continued) 
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Column No . Value/Format 

5 1 (15) 

10 1 ( 110 ) 

8 fields of 10 ElO.O 

10 2 ( 110 ) 

8 fields of 10 ElO.O 

10 3 ( 110 ) 

8 fields of 10 ElO.O 

10 4 (110) 

8 fields of 10 ElO.O 

10 5 (110) 

8 fields of 10 ElO.O 

10 6 ( 110 ) 

8 fields of 10 ElO.O 

10 7 (110) 

8 fields of 10 ElO.O 

10 8 ( 110 ) 

8 fields of 10 ElO.O 

10 9 (110) 

8 fields of 10 ElO.O 

1,2 10 ( 110 ) 

8 fields of 10 v ElO.O 

Figure 2. Input Names and Formats 


Parameter Name/Definition - 

NTS = number of bearing support 
points 

Indicates that the ball diameters 
follow 

D(I) (I = 1 to NTS) 

Indicates inner raceway diameter 
follows 

Did) 

Indicates outer raceway diameter 
follows 

DO (I) 

Indicates ball modulus of elasticity 
follows 

El (I) 

Indicates race modulus of elasticity 
follows 

E2 (I) 

Indicates shaft modulus of elasticity 
follows 

EH(I) 

Indicates bearing contact angle 
follows 

ALPHA (I) 

Indicates temperature coefficient 
for viscosity follows 

* ALPHAT(I) 

Indicates pressure coefficient for 
viscosity follows 

* ALPHMU (I) 

Indicates the density of the lubri- 
cant follows 

DENSE (I) 

Page 1 of 3 
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Column No . Value/Format 

9,10 11 (110) 

8 fields of 10 ElO.O 

9,10 12 (110) 

8 fields of 10 ElO.O 

9,10 13 (110) 

8 fields of 10 ElO.O 

9,10 14 (110) 

8 fields of 10 ElO.O 

9,10 15 (110) 

10 1 ( 110 ) 

8 fields of 10 ElO.O 

9,10 16 (110) 

8 fields of 10 ElO.O 

9,10 17 (110) 

8 fields of 10 ElO.O 

9,10 18 (110) 

8 fields of 10 ElO.O 

9,10 19 (110) 

8 fields of 10 ElO.O 

9,10 20 (110) 

8 fields of 10 ElO.O 

Figure 2 (Continued) 


Parameter Name/Definition 

Indicates lubricariit kinematic 
viscosity follows 

NUO(I) 

Indicates ball Poisson's ratio 
follows 

NUl(I) 

Indicates raceway Poisson's ratio 
follows 

NU2 (I) 

Indicates the bearing basic static 
capacity follows 

CS(I) 

Indicates the number and value of 
concentrated loads to the shaft 
follow 

Temporarily limited to "1" 

PH(I,K) (K = 1) 

Indicates the shaft moment of 
inertia follows 

IH(I) 

Indicates the ball density follows 
DENSE (I) 

Indicates the shaft span follows 
LH(I) 

Indicates lubricant thermal con- 
ductivity 

KH(I) 

Indicates the RPM of the rotating 
raceway follows 

NH(I) 

Page 2 of 3 
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Column No. 


Value/Format 
9,10 21 (110) 

8 fields of 10 ElO.O 

9,10 22 (110) 

10 

8 fields of 10 
9,10 

8 fields of 10 110 

9,10 24 (110) 

8 fields of 10 110 

8 fields of 10 ElO.O 

9,10 25 (110) 

10 110 

9,10 26 (110) 

10 1 (110) 

8 fields of 10 ElO.O 

Blank Card 

1 through 5 99999 


Parameter Name/Definition 

Indicates bearing thrust load 
follows 

FAH(I) 

Indicates the number of applied 
torques to the shaft and their 
values follow 

Temporarily limited to "1" 

TH(I,K) (K=l) 

Indicates bearing number of balls 
follows 

ZH(I) 

Indicates lubricant type and 
sliding friction coefficients follow 

LTYPE(I) (0,1, 2, 4 indicates none, 
mist, grease, jet, respectively) 

MU(I) 

Indicates the pointer for rotating 
raceway follows 

IRRC (0,1 for outer, inner, 
respectively) 

Indicates the number and value of 
shaft torque arm lengths follow 

Temporarily limited to "1" 

AH(I,K) (K=l) 

Indicates end of case 


Indicates End of Jpb 


* Indicates this data will be determined from Table 
interpolation at some later date. 


Figure 2 (Continued) Page 3 of 3 


1 ( 110 ) 

ElO.O 

23 (110) 
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APPENDIX 


COMPUTER LISTING OF BABHAP 
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